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A recent updated version of the Lanzhou quantum molecular dynamics is reviewed, in which the 
momentum dependence of symmetry potential and the effective mass splitting of proton and neutron 
in nuclear medium are included in the model. The in-medium nucleon-nucleon (NN) elastic cross 
sections are evaluated with the scaling approach according to the effective mass and their influence on 
collective flows in heavy-ion collisions are discussed. The inelastic cross sections, in particular for the 
process of the channel NA — >-NN, are parameterized in accordance with the available experimental 
data. It is found that the in-medium cross sections play a signiflcant role on isospin emissions and 
result in a flat distribution for transverse flows and elliptic flows of free nucleons comparing with the 
in- vacuum ones. The rapidity distribution of the difference between neutron and proton transverse 
flows is sensitive to the stiffness of nuclear symmetry energy as a promising observable, which can 
not be influenced by the in-medium effect and collision centrality. Furthermore, the elliptic flow of 
free energetic nucleons in the case of the mass splitting of rrin > rrip is also related to the symmetry 
energy. However, the pion flows weakly depend on the symmetry energy and the mass splitting. 
PACS number(s): 21.65.Ef, 24.10.Lx, 25.75.-q 



I. INTRODUCTION 



The high-density behavior of nuclear symmetry energy 
is of importance in understanding the properties of com- 
pact stars, in particular, the structure of neutron stars 
(mass and radius), the cooling of protoneutron stars, the 
nucleosynthesis during supernova explosion of massive 
stars etc (HQ. Up to now, the high-density informa- 
tion of the symmetry energy is poorly known ■ The 
importance of high-density symmetry energy in hadron- 
quark phase transition was stated and investigated thor- 
oughly by Di Toro et al. [6|]. Heavy- ion collisions in- 
duced by neutron-rich beams at intermediate energies 
are an unique approach for extracting the information 
of isospin-asymmetric nuclear equation of state (EOS) 
in terrestrial laboratories. The symmetry energy in 
neutron-rich matter is usually expressed through the en- 
ergy per nucleon as E{p, 5) = E{p, S = 0) + i?sym(p)<5^ + 
0{6^) in terms of baryon density p = p„ + pp, relative 
neutron excess 6 — (pn—Pp) / (Pn+Pp), energy per nucleon 
in a symmetric nuclear matter E{p, ^ = 0) and nuclear 
symmetry energy per nucleon i?sym- From that, one 
can easily get the physical quantities of nuclear matter, 
such as pressure, chemical potential and single particle 
potential etc. Based on several complementary analysis 
of available experimental data associated with transport 
models, a symmetry energy of Esym{p) ~ 31.6(/o/po)^ 
MeV with 7 = 0.69 — 1.05 was extracted for densities be- 
tween O.lpo and 1.2po iJf- However, predictions for the 
high-density symmetry energy based on various micro- 
scopical orphenomenological many-body theories diverge 
widely [7H9j. Moreover, the high-density behavior of nu- 
clear symmetry energy predicted by transport models as- 
sociated with the existing experimental data are also dif- 
ferent largely because of the inconsistent treatment of the 



mean-field potentials and the unclear in-medium nuclear 
interaction and properties of resonances and an 

opposite conclusion was drawn. Furthermore, new ex- 
perimental data related to high-density phase diagram 
and modifications of transport models including the in- 
medium effects in two-body collisions and in mean-field 
propagation, are very necessary. 

Based on the QMD ('quantum' molecular dynamics) 
approach over the past decade the LQMD (Lanzhou 
quantum molecular dynamics) model has been success- 
fully updated and applied to treat the nuclear dynamics 
from near Coulomb barrier energies up to relativistic en- 
ergies (several A GeV) [lil - fisj . in which the isospin ef- 
fect, nuclear structure (shell effect) in mean- field poten- 
tial, density dependence of symmetry energy, fermionic 
nature of nucleon in mean-field evolutions and in NN col- 
lisions, and inelastic channels including the in-medium 
effect etc, have been improved thoroughly. Furthermore, 
the momentum dependence of the symmetry potential 
was also included in the model, which results in an ef- 
fective mass splitting of proton and neutron in nuclear 
medium and its influence on isospin emission was investi- 
gated . In this work, we present a review of the recent 
updated LQMD model. Collective flows of free nucleons 
(light clusters) and pions produced in heavy-ion collisions 
and influence of the in-medium effect such as NN cross 
section, optical potential in nuclear medium etc, colli- 
sion centrality, mass splitting and stiffness of symmetry 
energy are investigated systematically. Extraction of the 
high-density symmetry energy from flow information are 
discussed. 

The article is organized as follows. In Sec. II we give 
a description of the recent version of LQMD model. In- 
medium modiflcation and its influence on isospin particle 
emissions, collective flows, and some promising observ- 
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ables for extracting the high-density symmetry energy 
are presented in Sec. III. In Sec. IV conclusions are 
discussed. 



II. MODEL DESCRIPTION 

It is well known that the QMD model has the advan- 
tage of correctly treating the A^-body correlations caused 
by the overlapping of the Gaussian wave packets and by 
the event- by- event simulations. A similar Gaussian wave 
packet has been used in the stochastic mean-field (SMF) 
model and obtained progress in the description of nuclear 
multi-fragmentation reactions by Catania group 0, [lO] ■ 
The same as in the QMD [l^l ; the wave function for each 
nucleon or resonance in the LQMD model is represented 
by a Gaussian wave packet. The width of the wave packet 
is fixed in the mean-field evolutions, but increases with 
the mass number of nuclide [l6l - [l8j . The total A^-body 
wave function is assumed as the direct product of the co- 
herent states and the anti-symmetrization is neglected. 
Therefore, the fermionic nature of nucleon is lost in the 
QMD-like models. We used a phase constraint approach 
to embody the nucleonic Fermi properties [l6| . From 
the Wigner transformation for the total wave function, 
one can get the density distributions in coordinate and 
momentum space. In the LQMD model, the time evolu- 
tions of the baryons (nucleons and resonances (A(1232), 
N*(1440), N*(1535))) and mesons in the system under 
the self-consistently generated mean-field are governed 
by Hamilton's equations of motion, which read as 



The local potential energy-density functional reads 
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(1) 



Here we only consider the Coulomb interaction for 
charged hyperons. The Hamiltonian of baryons consists 
of the relativistic energy, the effective interaction poten- 
tial and the momentum dependent interaction. The ef- 
fective interaction potential is composed of the Coulomb 
interaction and the local potential 

Uint — Ucoul + Uloc- (2) 

The Coulomb interaction potential is written as 

Ucoui = ^ E —erfinj/ViL) (3) 

where the is the charged number including protons 
and charged resonances. The r,j — jr^ — | is the relative 
distance of two charged particles. 

The local interaction potential is derived directly from 
the Skyrme energy-density functional and expressed as 
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where the p„, pp and p — Pn + Pp are the neutron, 
proton and total densities, respectively, and the S = 
{pn — Pp)/{Pn + Pp) is the isospin asymmetry. The coef- 
ficients a, j3, 7, gr, gsur, <?sf,° are related to the Skyrme 
parameters t(j,ti,t2,t^ and xq, Xi, X2, [l^ and the pa- 
rameter Sly6 is taken in the calculation with the values of 
9.9 MeV, 22.9 MeV fm^, -2.7 MeV fm^ and 0.16 frn'^ for 
<7r, gsur, gYur ^^'^ Po respectively. The bulk parameters 
a, P and 7 are readjusted after inclusion the momentum 
term in order to reproduce the compression modulus of 
symmetric nuclear matter (here, K=230 MeV) and the 
binding energy of isospin symmetric nuclear matter at 
saturation density, which have the values of -296.6 MeV, 
197 MeV and 1.143, respectively. The E[°yl^ is the local 
part of the symmetry energy, which can be adjusted to 
mimic predictions of the symmetry energy calculated by 
microscopical or phenomenological many-body theories 
and has two- type forms as follows: 
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Es°ym{p) = 7^Csym{p/ PoV% 



and 



KZnip) = O.sy,n{p/Po) + bsym{p/po) 



(6) 



(7) 



The parameters Csynn o,sym and hsym are taken as the 
values of 52.5 MeV, 43 MeV, -16.75 MeV and 23.52 MeV, 
32.41 MeV, -20.65 MeV corresponding to the mass split- 
tings of m* > m* and m* < m*, respectively. The 
values of 7^=0.5, 1., 2. have the soft, linear and hard 
symmetry energy with baryon density, respectively, and 
the Eq. (7) gives a supersoft symmetry energy, which 
cover the largely uncertain of nuclear symmetry energy, 
particularly at the supra-saturation densities. The local 
part of the symmetry energy can be adjusted to reflect 
the largely uncertain behavior of the symmetry energy at 
sub- and supra-normal densities. 

A Skyrme-type momentum-dependent potential is 
used in the model as follows: 



x/,(r,p,t)[ln(e(p- p')' + l)]2/,(r, p', t). (8) 



The term is also given from the energy- density functional 
in nuclear matter as 
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x/,(r,p)[ln(e(p - p')' + l)]Vr'(r, p')- (9) 
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Here Ct,t — CmomiX "I" ^t.t' — Cmom{^ ~ (''' 7^ ''"') 

and the isospin symbols r(r') represent proton or neu- 
tron. The sign of x determines different mass split- 
ting of proton and neutron in nuclear medium, e.g., 
positive signs corresponding to the case of m* < m*. 
The parameters Cmom and e was determined by fitting 
the real part of optical potential as a function of inci- 
dent energy from the proton-nucleus elastic scattering 
data. In the calculation, we take the values of 1.76 
MeV, 500 c^/GeV^ for the Cmom and e, respectively, 
which result in the effective mass m* /m=0.75 in nu- 
clear medium at saturation density for symmetric nu- 
clear matter. The parameter x can be adjusted as the 
strength of the mass splitting, and the values of -0.65 
and 0.65 are respective to the cases of m* > m* and 
TO* < m*. One can get the contribution of momentum- 
dependent interaction to symmetry energy for a cold nu- 
clear with Fermi distribution for the phase-space den- 
sity (r , p) = (r ) e (pF (r ) - 1 p I ) / (47rp3, (t) /3) with the 
Fermi momentum p;'(r) — h{3TT^ prY^^ 19]. The sym- 
metry energy per nucleon in the LQMD model is com- 
posed of three parts, namely the kinetic energy, the local 
part and the momentum dependence of the potential en- 
ergy as 
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After an expansion to second order around the normal 
density, the symmetry energy can be expressed as 
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Esym{p) ~ Esym{po) + — 
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slope parameter of L = 
and a curvature parameter of 
sym - dPo{9'^Esym/dp'^)\p=po. The values of slope 
parameters are 203.7 MeV, 124.9 MeV, 85.6 MeV 
and 74.7 MeV for the hard, linear, soft and supersoft 
symmetry energies, respectively, and the corresponding 
448 MeV, -24.5 MeV, -83.5 MeV and -326 MeV for 
the curvature parameters. Figure 1 is a comparison 
of different stiffness of nuclear symmetry energy after 
inclusion of the momentum-dependent interactions. All 
cases cross at saturation density with the value of 31.5 
MeV, which cover the largely uncertain of symmetry 
energy in a large density range. 

Combined Eq. (5) and Eq. (9), we get a density, 
isospin and momentum dependent single-nucleon poten- 
tial in nuclear matter as follows: 
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FIG. 1. (Color online) Density dependence of the nuclear 
symmetry energy for the cases of supersoft, soft, linear and 
hard trends with MDI interaction. 
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Here r ^ t' , dS'^/dpn = ASpp/p"^ and dS'^/dpp = 
—ASpn/p^- The nucleon effective (Landau) mass in nu- 
clear matter of isospin asymmetry 5 ~ (pn ^ Pp)/ {pn + Pp) 
with Pn and pp being the neutron and proton den- 
sity, respectively, is calculated through the potential as 
TO* = TUr/ (l + jpj-|^l) with the free mass to^. at 
Fermi momentum p = p^,. Therefore, the nucleon ef- 
fective mass only depends on the momentum-dependent 
term of the nucleon optical potential. The isovector 
part of the optical potential, i.e., the symmetry or Lane 
potential 2l|, can be evaluated from the expression 

U sym 

(p,p) = (t/„(p, 5, p) — C/p(p, (5, p))/2^. Shown in 
Fig. 2 is a comparison of the momentum dependence 
of single-nucleon optical potential with the mass split- 
tings of TO* > TO* in the left window and to* < to* in 
the right window for the hard and supersoft symmetry 
energies, respectively. One should note a cross appears 
in the case of to* > to* at high momentum, but which 
does not take place in the mass splitting of to* < to* and 
a broader separation exists with increasing the nucleon 
momentum. The symmetry potential is also being influ- 
enced by the mass splitting and just an opposite trend 
appears as shown in Fig. 3. A pronounced difference 
of both mass splittings can be constrained from elliptic 
flow in heavy-ion collisions, especially from the momen- 
tum (kinetic energy) distribution of Vj" ~ ^2 mid- 
rapidity [i^]. More specifically, the symmetry potential 
is plotted as functions of momentum and density in Fig. 
4. Both the stiffness of symmetry energy and the mass 
splitting affect the spectrum. The effective mass is oppo- 
site in both mass splittings as shown in Fig. 5, which is 
pronounced in the \A GeV energies (l~3po). Therefore, 



4 



both of the stiffness of nuclear symmetry energy and the 
mass splitting affect the reaction dynamics, in particular 
for the isospin emission in heavy-ion collisions. 

The same with the QMD or BUU (Boltzmann-Uehling- 
Uhlenbeck) model, a hard core scattering in two-particle 
collisions is assumed by using Monte Carlo procedures, in 
which the scattering of two particles is determined by a 
geometrical minimum distance criterion d < y^O.latot/''^ 
fm weighted by the Pauli blocking of the final states 
[3, IS]- Here, the total NN cross section atot in mb 
is the sum of the elastic and all inelastic cross sec- 
tions. The probability reaching a channel in a collision 
is calculated by its contribution of the channel cross sec- 
tion to the total cross section as Pch = (^ch/crtot- The 
choice of the channel is done randomly by the weight 
of the probability. We parametrized the total, elastic 
and inelastic NN cross sections in accordance with the 
available experimental data in free space as shown in 
Fig. 6. One expects the in-medium NN elastic col- 
lisions are reduced in comparison with the free-space 
ones. The in-medium elastic cross section is scaled ac- 
cording to the effective mass that was used in the BUU 
model [25|]. Following the IBUU04 model [111, the elastic 
cross section in the nuclear medium is evaluated through 

^■medium ^ {n^, ^ / nY CT^N ^^1*^ ^'^^ f^*NN and fiNN 

being the reduced masses of colliding nucleon pairs in 
the medium and in the free space, respectively. Shown in 
Fig. 7 is a comparison of the scaling factor as functions 
of baryon density, nucleon momentum and isospin asym- 
metry for the different mass splittings of to* > to* and 
TO* < TO*, respectively. It is interesting to notice that 
a splitting of nn (neutron- neutron), pn (proton- neutron) 
and pp (proton-proton) exists because of the mixing of 
neutron and proton effective masses. The mass splitting 
results in an opposite splitting of the in-medium cross 
sections of nn and pp channels. Recent microscopic cal- 
culations from the Brueckner-Hartree-Fock and the rela- 
tivistic Brueckner-Hartree-Fock basically favor the varia- 
tion of the in-medium cross section with baryon density, 
nucleon momentum and isospin asymmetry |27l [28| . 

The primary products in nucleon-nucleon (NN) colli- 
sions in the region of 1^ GeV energies are the resonances 
A(1232), Af*(1440), A^*(1535) and the pions. We have 
included the reaction channels as follows: 

NN^NA, NN^NN*, NN ^ AA, 
A o Nn, N* ^ Nn, NN NNn{s - state), 
A^*(1535) Nt]. (13) 



At the considered energies, there are mostly A reso- 
nances which disintegrate into a n and a nucleon in 
the evolutions. However, the N* yet gives considerable 
contribution to the energetic pion yields. The energy 
and momentum-dependent decay widths are used in the 
model for the A(1232) and iV*(1440) resonances [III- 
We have taken a constant width of r=150 MeV for the 



TV* (1535) decay. The in-medium effect of inelastic cross 
section is of importance in the estimation of meson pro- 
duction, which has been an interesting subject and still 
is an open p roblem because of the coupling of each chan- 
nels |29l43l| . A parameterized cross section is used in 
the LQMD model for the channel of NA -J>NN by fit- 
ting the available experimental data 3^ with a formula 
by Wolf et at [33], and compared with the cross section 
from the NN— >NA channel using the detailed balancing 
principle as shown in Fig. 8. The in-medium modifica- 
tions of the absorption of A in heavy-ion collisions are 
also investigated in Refs. H^. Calculations based 
on relativistic Dirac-Brueckner also favor the decrease 
trend of the A absorption cross section with increasing 
energy in nuclear medium at low densities (less than 2po) 
(36| . We used the parametrized cross sections calculated 
by the one-boson exchange model [37| for resonance pro- 
duction and the absorption of N* with the detailed bal- 
ancing principle. The strangeness is created by inelastic 
hadron-hadron collisions. We included the channels as 
follows: 

BB BYK, BB BBKK, Bn YK, Bn NKK, 
Yn ^ BK, BK ^ Yn, YN KNN. (14) 

Here the B strands for (N, A, N*) and Y(A, E), K(K0, 
K+) and K{K^, K~). The elastic scattering between 
strangeness and baryons are considered through the 
channels KB ^ KB, YB ^ YB and KB -^KB. At 
the moment, we do not implement the charge-exchange 
reactions between the KN KN and YN — YN chan- 
nels, such as K°p — > K+n, K+n — > K^p, An — > S^p, 
Art — >■ E'^n etc, which would be important in evaluating 
the isospin ratios of K^ / K^ and for probing the 

high-density symmetry energy. 

Analogously to baryons, the evolution of mesons (here 
mainly pions and kaons) is also determined by the Hamil- 
tonian, which is given by 



Nm 



Here the Coulomb interaction is given by 



Coul 



Nb 



(15) 



(16) 



where the Nm and Nb are the total numbers of mesons 
and baryons including charged resonances. The pion op- 
tical potential Re^,^^* originates from the medium effects 
in the hot and dense nuclear matter. In the calcula- 
tion, we can also choose the value in vacuum, that is, 
the ReT^^^* is set equal zero. The influence of the pio- 
nic mean field in heavy-ion collisions on the transverse 
momentum distribution was investigated by using a phe- 
nomenological ansatz and a microscomc approach based 
on the A-hole model by Fuchs et al. 38| . Here we use the 
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FIG. 2. (Color online) Momentum dependence of single-nucleon optical potential for isospin symmetric matter (5=0) and 
neutron-rich matter (3=0.2) with the mass splittings of m* > m* (left panel) and m* < m* (right panel), respectively. 




FIG. 3. (Color online) Symmetry potential as a function of momentum with the hard and supersoft symmetry energies and 
with different mass splittings at nuclear matter density 0.5po, po and 2po, respectively. 



Dhenomenological ansatz suggested by Gale and Kapusta 
[39j. Then the dispersion relation reads 

^{P^, P^) = \/i\pJ^P0)^^^0 - U, (17) 

U = \/pI + 1^1- m^, (18) 

mo = TO,i- + 6.5(1 — a;^")m7r, (19) 
pI = {1- x)^ml + 2mom^{l- x). (20) 



The phenomenological medium dependence on the 
baryon density is introduced via the coefficient x{pi) = 
exp(— a(/9i//0o)) with the parameter a = 0.154 and the 
saturation density po in nuclear matter. Influence of the 
in-medium effects on the charged pion ratio is also inves- 
tigated in Ref. |4(|. We consider two scenarios for kaon 
(antikaon) propagation in nuclear medium, one with and 
one without medium modification. From the chiral La- 
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FIG. 4. (Color online) Contour of symmetry potential versus momentum and density of nuclear matter with different symmetry 
energies and mass splittings. 
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FIG. 5. Effective mass of proton and neutron in nuclear medium as a function of baryon density with the isospin asymmetry 
of 5=0.2 for the mass splittings of m* > m* (left panel) and m* < m* (right panel), respectively. 



grangian the kaon and antikaon energy in the nuclear 
medium can be written as 41, 43] 



1/2 



(21) 



w(Pi,Pi) 



'^^ + v1 -a-j^pf + {bKPi) 



1/2 



^KPi-, 

(22) 

respectively. Here the bx = 3/(8/^) wO.32 GeVfm^, 
the ttK and a-j^ are 0.18 GeV^fm^ and 0.3 GeV^fm^, re- 
spectively, which result in the strength of repulsive kaon- 
nucleon potential and of attractive antikaon-nucleon po- 



tential with the values of 25.5 MeV and -96.8 MeV at 
saturation baryon density. Shown in Fig. 9 is comparison 
of the pion self-energy and optical potential as a function 
of momentum. One notices that an attractive potential 
is used in the LQMD model, in particular for energetic 
pions in nuclear medium. We calculated the pion emis- 
sions after inclusion of all, part and none of the optical 
potential and Coulomb interaction of pions with charged 
baryons in the reaction -'^^^Sn-^-'^^'^Sn for near-central col- 
lisions (b=l fm) at incident energy of E;afc=400yl MeV 
as shown in Fig. 10. Without the Coulomb and pion 
potential, a flat tt^/tt^ ratio appears with the transverse 
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FIG. 6. Nucleon-nucleon cross sections of total, elastic and inelastic channels in free space parameterized in the LQMD model 
as a function of nucleon incident energy in the laboratory frame. 
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FIG. 7. In-medium nucleon-nucleon elastic cross section scaled by nucleon effective mass as functions of baryon density {S = 0.2), 
nucleon momentum (5 — 0.2) and isospin asymmetry (p — po). 



momentum. The Coulomb force enhances the energetic 
7r+ production because of its repulsive force with charged 
baryons, which is opposite for n~ emission, consequently 
increases the low-energy tt^/tt^ ratio and reduces the 
high-energy yields. However, the pion optical potential 
is always attractive for both charged pious. Calculations 
in Ref. [l^l show that the pion potential slightly affects 
the total tt^/tt^ yields owing to without distinguishing 
isospin effects for the charged pion potentials. Experi- 
mental measurements on the distributions of the trans- 



verse momentum or kinetic energy for the '!t~ /t:~^ ratio 
will be helpful in understanding the pion optical potential 
in nuclear medium. 

III. RESULTS AND DISCUSSIONS 

Collective flows in heavy-ion collisions have been stud- 
ied well to extract the information of EOS at extreme 
conditions, i.e., high density, high temperature and large 



8 



400 



300 



200 



100 




FIG. 8. Comparison of the cross section of the channel 
NA — j-NN parameterized in the model and compared with 
the approach by Bertsch and Gupta used in the BUU model 
and with the experimental data [32 |. 



isospin asymmetry etc., and to explore the knowledge 
of in-medium properties and nuclear dynamics, such as 
nucleonic flow, light cluster flow, and meson flow etc 

Ha 



[131 1431 I44j . It is possible to reconstruct the reaction 
plane with flow analysis and hence to study azimuthal 
correlations of emitted particles. The flow information 
can be expressed as the first and second coefficients 
from the Fourier expansion of the azimuthal distribution 
^{y^Pt) = iVo(l + 2Vi{y,pt) cos(0) + 2^2 cos(20)) 



|45| , where the pt — yPx+Py ^^'^ V ^'^^ the trans- 
verse momentum and the longitudinal rapidity along the 
beam direction, respectively. The directed (transverse) 
flow is defined as the first coefficient and expressed as 
Vi = {Px/Pt)i which provides the information of the az- 
imuthal anisotropy of the transverse emission. The ellip- 
tic flow V2 = {{px/ptY ~ {Py/PtY) gives the competition 
between the in-plane (V2 >0) and out-of-plane (V2 <0, 
squeeze out) emissions. The brackets indicate averaging 
over all events in accordance with a specific class such 
as rapidity or transverse momentum cut. Shown in Fig. 
11 is a comparison of rapidity distributions of transverse 
and elliptic flows of free nucleons for near-central (b=l 
fm) and semi-central (b=6 fm) ^^^Au-|-^^^Au collisions 
at the incident energy of 400^4 MeV with a choice of 
the mass splitting of m* > m* for the hard and super- 
soft symmetry energies, respectively. One notices that 
the difference of neutron and proton transverse flows is 
small for the near-central collisions, but pronounced for 
the semi-central collisions, in particular around projectile 
rapidity. The elliptic flow disappears for the near-central 
heavy-ion collisions and the influence of symmetry en- 
ergy on the elliptic flow distribution is very weak. The 
difference of neutron and proton elliptic flows is obvious 
in the semi-central collisions. In this work, the free nu- 



cleons and fragments are constructed with a coalescence 
model at freeze-out, in which nucleons of the reaction sys- 
tem are considered to belong to one cluster in the phase 
space with the relative momentum smaller than Pq E^nd 
with the relative distance less than Rq (here, Po=200 
MeV/c and Rq—?> fm). The set of the parameters can 
reproduce the distribution of experimental charged frag- 
ments in a large region of energies [1^ . 

It has been well known that the transverse flow sensi- 
tively depends on the variation of in-medium NN cross 
sections, especially around the balance energies 46. 47|. 
As a comparison, we calculated the transverse flow of free 
nucleons as shown in Fig. 12 for the "'^^^Sn-|--'^^^Sn reac- 
tion at the incident energy of 400 MeV/nucleon using the 
elastic cross sections in free space and in nuclear medium, 
respectively. One notes that the in-medium cross section 
leads to a flat distribution for the near-central collisions. 
However, its influence on the semi-central collisions is 
slightly. The in-medium cross section also has a strongly 
influence on the elliptic flow distribution as shown in Fig. 
13, which reduces the out-of-plane emissions of free nucle- 
ons, especially in the domain of mid-rapidity. Shown in 
Fig. 14 is the in-medium effect on light clusters with the 
hard (left window) and supersoft (right window) sym- 
metry energies. It is seen that the reduced in-medium 
cross sections slightly modify the collective flows of light 
clusters (^H and ^He). 

In order to extract the density dependence of symme- 
try energy, especially for the high-density information, 
we calculated the flow difference of free nucleons as a 
function longitudinal rapidity as shown in Fig. 15. One 
can see that the spectrum is sensitive to the stiffness of 
symmetry energy in both near-central and semi-central 
collisions, especially around the projectile rapidity. The 
NN cross sections in free space and in nuclear medium 
basically does not change the distributions. A supersoft 
symmetry energy has a more pronounced spectrum. As 
a comparison, we also present the transverse flow dif- 
ference of ■^H and ^He as shown in Fig. 16. A similar 
structure appears, but weakly depends on the symmetry 
energy. A narrow rapidity distribution is observed for the 
light cluster emission because a number of free nucleons 
already take away much kinetic energies dissipated from 
the reaction zone of colliding partners. Shown in Fig. 17 
is a comparison of transverse momentum distribution of 
elliptic flow of free nucleons with the hard and supersoft 
symmetry energies for the mass splittings of m* > m* in 
the left panel and of m* < m* in the right panel, respec- 
tively. The influence of symmetry energy on the elliptic 
flow is obvious only in the mass splitting of m* > m*, 
in particular at high transverse momentum. It should be 
mentioned that the distribution of elliptic flow is sensi- 
tively dependent on the mass splitting Constraint 
of the high-density symmetry energy from elliptic flow 
has been performed in experiment 13] and more precise 
measurements are planned at RIKEN and GSI. 
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FIG. 9. Momentum dependence of pion self-energy and optical potential at the baryon densities of 1.5po and 3po, respectively. 
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FIG. 10. Transverse momentum distributions of the ratios of 
7r~ to ■k'^ yields with including all, only Coulomb potential 
and none of the optical potential and Coulomb force of pions 
with baryons in the reaction ^■^'*Sn-|-^^*Sn at incident energy 
of E,,i,=400A MeV. 



Production of pions in heavy-ion collisions at near- 
threshold energies has been verified as a sensitive probe 
of nuclear symmetry energy, in particular the 7r~/7r+ ex- 
citation functions l3"12|, which can be easily detected 
in experimentally and be mainly created in the domain 
at supra-saturation densities of compressed nuclear mat- 
ter. Shown in Fig. 18 is a comparison of transverse flows 
of TT^ and 7r+ in ^^^Au and ^^""Au collisions at a near- 
threshold energy of 400 MeV/nuclcon with the hard and 
supersoft symmetry energies for different mass splittings 
of m* > m* in the left window and of m* < m* in the 
right window, respectively. It can be seen that the sym- 
metry energy and mass splitting basically does not affect 



the flow structure. Furthermore, an antiflow of 7r+ (solid 
and blank squares) is observed in the both mass splittings 
by comparing with the transverse flow of free nucleons in 
Fig. 11. The phenomena is caused by the fact that pions 
produced in the 'fire ball' are absorbed again by partic- 
ipant nucleons (shadowing effect). A more pronounced 
fiow (tt" ) or antiflow (tt^) is observed comparing with the 
collisions at a higher incident energy of 1.5 GeV/nucleon 
[l8| . The antiflow of 7r+ in heavy- ion collisions was also 
investigated using BUU Q and IQMD model For 
a comparison, we also calculated the flow difference be- 
tween 7r+ and tt^ with different symmetry energy and 
mass splitting as shown in Fig. 19. Although a pro- 
nounced flow difference {V^ — ) appears around the 
projectile rapidity, the influence of symmetry energy and 
mass splitting on the spectrum is still very weak. 



IV. CONCLUSIONS 

The influence of the in-medium effects in nuclear re- 
actions is investigated within the LQMD model, i.e., 
the NN cross sections, propagation of pions in nuclear 
medium. The in-medium cross sections are of impor- 
tance on isospin emissions and lead to a flat distribution 
on the transverse flows and elliptic flows of free nucleons 
comparing with the ones in free space. The rapidity dis- 
tribution of the transverse flow difference between neu- 
tron and proton depends on the density dependence of 
nuclear symmetry energy, in particular around the pro- 
jectile rapidity. The in-medium effect and collision cen- 
trality slight affect the flow difference. The elliptic flow 
of free nucleons is also sensitive to the symmetry energy, 
especially at the high momentum in the case of the mass 
splitting of TO* > TO*. The pion optical potential plays 
a significant role on the transverse momentum distribu- 
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asy-hard asy-supersoft 




FIG. 11. (Color online) Rapidity distributions ol directed and elliptic flows of free nucleons in the ^''^Au+^^^Au reaction at the 
incident energy of 400 MeV/nucleon for the near central (b=l fm) and semi-central (b=6 fm) collisions using the NN elastic 
cross section in free space with the mass splitting of m* > m* for the hard (left panel) and supersoft (right panel) symmetry 
energies, respectively. 
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FIG. 12. (Color online) Comparison of the transverse flow of free nucleons with in-vacuum and in-medium cross sections in 
the ^^^Sn-f ^^''Sn reaction at the incident energy of 400 MeV/nucleon for the near-central (b=l fm) and semi-central (b=6 
fm) collisions with the mass splitting of m* > m* for the hard (left panel) and supersoft (right panel) symmetry energies, 
respectively. 
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FIG. 13. (Color online) The same as in Fig. 12, but for the elliptic flow for semi-central collisions (b=5 fm). 
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FIG. 14. (Color online) Rapidity distributions of directed and elliptic flows of light clusters in the ^^^Sn-f'^^'^Sn reaction at a 
beam energy of 400^4 MeV for the semi-central (b=5 fm) collisions with the mass splitting of m* > m* for the hard (left panel) 
and supersoft (right panel) symmetry energies, respectively. 



tion of tt^/tt^ ratio. The transverse flow of charged pions 
weakly depend on the symmetry energy and also the mass 
splitting. An antiflow for 7r+ emission is pronounced in 
both mass splittings 
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